Resonance Patterns in a Stadium-shaped Microcavity 
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We investigate resonance patterns in a stadium-shaped microcavity around UckR ^ 10, where 
Uc is the refractive index, k the vacuum wavenumber, and R the radius of the circular part of the 
cavity. We find that the patterns of high Q resonances can be classified, even though the classical 
dynamics of the stadium system is chaotic. The patterns of the high Q resonances are consistent 
with the ray dynamical consideration, and appears as the stationary lasing modes with low pumping 
rate in the nonlinear dynamical model. All resonance patterns are presented in a finite range of kR. 
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I. INTRODUCTION 

As a central concept of the promising photonics tech- 
nology, recently microcavity has attracted much atten- 
tion due to its necessity for the generation and switch- 
ing of coherent light In a cylindrical and spherical 
dielectric cavities, light trapped by total internal reflec- 
tion can make very high Q modes, so-called whispering 
gallery modes (WGM). Lasing of WGM has been stud- 
ied extensively 0, Q 1^- Imposing small deformation 
on circular and spherical symmetric cavities, one can ob- 
tain lasing modes with directional emission which is more 
useful for the potential application to photonics 0, 0- 

It is known that even in chaotic microcavities the las- 
ing modes with high Q and very good directionality can 
be excited 1^]. There are several reports for observation 
of scarred lasing modes in various microcavities whose 
patterns show enhanced field intensity alongan unsta- 
ble periodic orbit of the cavity |3 [TE El El. More- 
over, due to the inherent properties of the dielectric cav- 
ities, the existence of quasi-scarred resonance modes [T^j. 
which are not supported by any unstable periodic or- 
bit, are suggested and numerically confirmed in a spiral- 
shaped dielectric microdisk. In a nonlinear dynamical 
model of a stadium-shaped cavity with an active medium, 
Harayama et al. examined laser action on a single spa- 
tially chaotic wave function 14] and locking of two reso- 
nance modes of different symmetry classes and slightly 
different frequencies 15j. 

In this paper we study resonances of the stadium- 
shaped microcavity by using the boundary element 
method (BEM) 16]. We present the whole resonances 
in the range of 8 < Iie{nckR) < 11 with the refractive 
index nc = 2 to show what resonance patterns are pos- 
sible to appear as lasing modes and to find the relation- 
ships between Q factor and resonance pattern. We note 
that the classical dynamics of the stadium-shaped bil- 
liard is chaotic. This means that the resonances cannot 



be classified by mode indices, similar to the absence of 
good quantum number in quantum mechanical problem. 
In spite of its chaotic ray dynamical properties, we find 
that the high Q modes in this Re{nckR) range have spe- 
cial patterns which can be classified from other modes, 
and these modes are really excited with low pum ping rate 
in the nonlinear dynamical model developed in |l^ . 

This paper is organized as follows. In Sec ^ we sum- 
marize the BEM used to obtain resonances numerically, 
and the resulting resonances are presented and discussed 
in Sec lIIII In Sec lIVI the variations of the resonance po- 
sition and pattern with the refractive index ric are dis- 
cussed. We confirm, in Sec0 the fact that the high 
Q resonances are the stationary modes with low pump- 
ing rate in the nonlinear dynamical model, and a ray 
dynamical consideration for the high Q resonance pat- 
terns is given in Sec I VII Finally we summarize results in 
SeclVTTI 



II. BOUNDARY ELEMENT METHOD 

In this section we briefly describe the boundary el- 
ement method (BEM) 16] which is effective when the 
boundary is strongly deformed from a circular shape. In 
the microdisk case. Maxwell's equations can be reduced 
to a two-dimensional Helmholtz equation 



(V'+n(r)2A:2)V; = 0. 



(1) 
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The main idea of BEM is to replace the Helmholtz equa- 
tion by one-dimensional boundary equations and then to 
discretize the boundaries. To obtain the one-dimensional 
boundary equation, the Green function, the solution of 
the Helmholtz equation with a delta function source, and 
Green's identity are used. The Green function for the 
Helmholtz equation is given by the zeroth-order Hankel 
function of the first kind, 

G{r,r';nk) = -'-Hi'\nk\r-r'\). (2) 

Then the resulting boundary integral representation is 

cV'(r) = ^ ds' [^^,G{r, s'; nk)ip{s') - G(r, s'; nk)d^i}{s')], 

(3) 




-0.145 



-0.15- 




FIG. 1: Convergence of the solution kR for the resonance 
OO3 with the number of elements on the boundary, (a) Real 
part (b) Imaginary part 



where c is 1 for r inside the boundary F, 1/2 for r on the 
boundary, otherwise, and 9^ means the outward nor- 
mal derivative at the boundary. In our case the bound- 
ary is a simple stadium shape, we can therefore make, 
by locating r on the boundary, two boundary equations 
corresponding to the inside and outside of the microdisk, 
respectively. Discretization of the two boundary equa- 
tions leads to a matrix equation. 



M 



0, 



(4) 



where (ip) and (dui^) are column vectors composed by the 
values of i/j and dui^ at the element points on the bound- 
ary. The non-trivial solution of this matrix equation ex- 
ists only when the determinant of M is zero, which is the 
condition used in determining the complex wavenumber 
k of resonances. 

Note that Hankel functions in the Green function and 
its normal derivative are singular at the origin, so a 
careful treatment is needed to overcome this sigularities. 
The sigularity of h[^\z) in d^^G can be compensated 
by a geometrical factor, and the corresponding diagonal 
terms can be expressed in terms of the curvature hi of 
the boundary. In order to treat the sigularity of H^\z) 

in G, we use the asymptotic expression of H^\z) and 
integrate it in the related elements 16]. 

In this paper, we focus on TM polarization where both 
the wavefunction ip and its normal derivative dyij) are 
continuous across the boundary, and the ijj corresponds 
to z-component of electric field Ez when the disk lies on 
x-y plane. In the case of TE polariztion, the represents 
the magnetic field Hz and the wavefunction ip is continu- 
ous across the boundary, but its normal derivative is not, 
instead, n{r)~'^ dj^tjj is continuous. This boundary condi- 
tion can be incorporated into the above formalism by a 
slight correction. In practical calculations we take into 
consideration of the symmetry of the stadium-shaped mi- 
crocavity, i.e., stadium has symmetries about the x- and 
7/-axes. The symmetry consideration leads to the reduc- 
tion of computing time by allowing to treat only quarter 
part of the stadium boundary. Instead we have to use 
the desymmetrized Green function when constructing the 
matrix equation Eq.(^. 
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FIG. 2: (Golor online) Resonance points in the complex kR 
space. Black, red, blue, and green points represent resonances 
of odd-odd, even-odd, odd-even, and even-even symmetric 
class, respectively. The numbering starts from the smallest 
Ke{kR) value. 



Generally the space between element points can be 
taken as a variable of the curvature, i.e., if the curva- 
ture is large, then take a large number of points. The 
stadium-shaped boundary is, however, simple enough 
to take equal spacing points. In the practical calcu- 
lation we take 50 elements on the quarter part of the 
boundary, which corresponds to about 12 element per a 
wavelength(Ac = 27r/nck). Fig.l shows the convergence 
of the complex k value of the resonance OO3 (see Fig. 2) 
with the number of elements N in the quarter part of the 
boundary. The values at = 50 gives an error less than 
0.2 % from the limit values. 



III. RESONANCES IN STADIUM-SHAPED 
MICROCAVITY 

We find 33 resonances in the range 8 < UckR < 11 
with Tic = 2. The number of resonances in this range can 
be inferred from the modified Weyl's theorem jl 71]. 



k{s) 



(5) 



where A is the area, L the length of the boundary, and the 
-(+) refers to Dirichlet (Neumann) boundary conditions. 
If we take the Dirichlet boundary condition the num- 
ber of eigenmodes over the concerned range is given by 
A^o(ll) — ^0(8) — 30.0 and for Neumann boundary con- 
dition 7Vo(ll) - No{S) 0^ 34.78. Since under the bound- 
ary condition of the present problem both and du'p 
on the boundary are nonzero, it is resonable to find the 
number of resonances between those predicted in both 
boundary conditions. For each symmetry class we find 
7 (odd-odd), 8 (odd-even), 8 (even-odd), and 10 (even- 
even) resonances, respectively. The smaller number of 
resonances for the odd-odd symmetry case and the larger 
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FIG. 3: (Color online) Resonances with odd-odd symmetry, 
(a) OO4, kR = (4.67, -0.048) (b) OO7, /ci? = (5.35, -0.042) 



(c) OO2, (4.31, 
(e) OO3, A;i?= (4.49, 
(g) OOi, (4.07, 



-0.086) (d) OO5, kR = (4.91, -0.076) 
-0.14) (f) OOe, /ci? = (5.13, -0.092) 
-0.24) 



FIG. 4: (Color online) Resonances with even-odd symmetry, 
(a) EO3, kR = (4.33, -0.053) (b) EOq, kR = (5.01, -0.044) 
(c) EOi, kR = (4.01, -0.075) (d) EO4, kR = (4.61, -0.084) 
(e) EO7, kR = (5.23, -0.073) (f) EO5, kR = (4.80, -0.11) 
(g) EOs, kR = (5.46, -0.083) (h) EO2, kR = (4.25, -0.20) 



number for the even-even symmetry case, can be under- 
stood by Eq.(|5|), since the former is close to the Dirichlet 
boundary condition and the latter is close to the Neu- 
mann boundary condition. 

The relation between Q factor and the complex 
wavenumber k is given by 



Q = 



Re[kR] 
'2Im[kRy 



(6) 



Therefore, the smaller absolute value of imaginary part 
of k means the higher Q modes. Fig. 2 shows all reso- 
nance positions in the concerned kR range. Black points 
are the resonances with odd-odd parities on x and y co- 
ordinates, respectively. This symmetry class is denoted 
by OOi{i = 1,2,---) and the index numbering starts, 
for convenience, from the resonance with the smallest 
Ke{kR) in the concerned range. With the same way 
we label other resonances, and in the figure with differ- 
ent color we distinguish the different symmetry classes, 
red, blue, and green points denote even-odd, odd-even, 
and even-even parities, respectively. Unlike the circular 
boundary in which WG modes have very small absolute 
value of lm{kR) ( e.g., -0.0035 < lm{kR) < -0.0025 
on Re(/ci?) 8), in the stadium-shaped microcavity 
the high Q modes have relatively large absolute value. 



\lm{kR)\ :^ 0.04. Remember that the classical dynam- 
ics of the stadium billiard is chaotic, so it is impossible 
for ray to go around infinitly trapped by total internal 
reflection like WGM. This fact explains the absence of 
resonances with \lm{kR)\ < 0.03. 

It is interesting to see the relationship between pat- 
terns and lm{kR) values of resonances. The resonance 
patterns corresponding to all resonance points in Fig. 2 
are drawn in Fig. 3 (odd-odd parity). Fig. 4 (even-odd 
parity). Fig. 5 (odd-even parity), and Fig. 6 (even-even 
parity). In figures color changes in order of dark brue- 
brue-white- red-dark red as the intensity of electric field 
increases. From Figs. 3-6, several resonances have similar 
patterns. In particular the parity on y coordinate plays 
more crucial role in classifying the patterns rather than 
the parity on x coordinate. 

First let us focus on the high Q resonances for all 
symmetry classes. OO4, 007(Fig.3 (a),(b)) and EO3, 
£^06(Fig.4 (a),(b)) are high Q modes with odd-odd and 
even-odd parities, respectively. It is obvious that these 
resonances can be classified as one class of pattern by 
observing the distribution of high intensity field near the 
boundary and same directionality of near field emission. 
The only difference between them is the number of in- 
tensity spots along the boundary. This fact can be il- 



4 



(a) (b) 




FIG. 5: (Color online) Resonances with odd-even symmetry, 
(a) OEi, kR = (4.25, -0.040) (b) OE5, kR = (4.89, -0.055) 
(c) OE2, kR = (4.37, -0.094) (d) OEq, kR = (5.05, -0.075) 
(e) OE3, kR = (4.58, -0.10) (f) OEg, kR = (5.32, -0.096) 
(g) OE7, kR = (5.20, -0.18) (h) OE4, kR = (4.83, -0.25) 




(a) (b) 




lustrated by almost constant difference of Re(/ci?), i.e., 
ARe{kR) between adjacent two modes among these is 
about 0.34. They show very good directionality of emis- 
sion. Similar discussion can be applied to the high Q 
modes with odd-even and even-even parities. OEi^ OE^ 
(Fig.5 (a),(b)) and EE3, EEs (Fig.6 (a),(b)) are high 
Q modes with odd-even and even-even parities, respec- 
tively. From the fact that the difference of Ke{kR) be- 
tween succesive modes are ARe{kR) o:^ 0.31, we can con- 
firm that these modes can be classified as a class. The 
resonance pattern of this class is somewhat different from 
the high Q modes of odd parity on y coordinate, high in- 
tensity spots are shown on the semicircular part of the 
boundary and on the central part of the cavity. We note 
that these patterns have high intensity on the central 
part. In a circular microdisk the high Q modes are WG 
modes which show high intensity near the boundary, indi- 
cating that the WG modes are supported by rays circulat- 
ing and trapped by the total internal reflection. Although 
the present stadium-shaped microdisk shows chaotic ray 
dynamics, implying no exact WG modes trapped by total 
internal reflection, similar patterns are shown in the for- 
mer high Q modes such as OO4, OO7, £^03, EOq. There- 
fore, the high intensity on the central part of the cavity in 
the latter high Q modes {OEi, OE^, EEs, EEs) means 



FIG. 6: (Color online) Resonances with even-even symmetry, 
(a) EEs, kR = (4.57, -0.044) (b) EEs, kR = (5.20, -0.065) 
(c) EEi, kR = (4.04, -0.11) (d) EE4, kR = (4.71, -0.083) 
(e) EEq, kR = (5.38, -0.069) (f) EE2, kR = (4.23, -0.10) 
(g) EEg, kR = (4.95, -0.10) (h) EEio, kR = (5.47, -0.13) 
(i) EEr, kR = (4.98, -0.22) (j) EE5, kR = (4.74, -0.27) 



that rays supporting those modes are different from those 
for the WG type resonances, and they would pass the cen- 
tral part of the cavity. The ray dynamical consideration 
for the high Q modes would be presented in Sec lVll 

Next we consider other pattern classes of resonances 
in the same Re{kR) range. EOi, OO2, EO4, 00^, and 
£^07 form a class of pattern which has attached two cir- 
cle structure inside the boundary and ARe(/ci?) 0.3. 
EEi^ OE2, EE^^ OEq^ and EEg show similar two circle 
structure with one line of spots near the straight seg- 
ments of the boundary, giving ARe{kR) :^ 0.33. EE2, 
OEs, EEq, and OEs have a simpe line structure, the 
high intensity spots are concentrated on the x-axis, giv- 
ing ARe{kR) 0.36. OO3, EO3, OOq, and EOs 
show low intensity on the center of the stadium, giving 
ARe{kR) ^ 0.32. The very high loss (large \lm{kR)\) 
resonances are OOi, £^02, EEr^^ OE^^ EEj and OEj 
which are basically bouncing ball type resonances. 
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FIG. 7: (Color online) Variation of resonance with ric. (a) 
Variation of Re{kR) of the resonance mode OO4. (b) Varia- 
tion of Im(kR) of the resonance mode 004- The inset show 
the convergence of Im(A;i?) to a nonzero value clearly when 
iV = 50. (c) The limiting resonance of the mode OO4 (at 
no/ric = 0.01). (d) The limiting resonance of the mode OO5 
(at no/ric = 0.01). 



From the patterns we can find general rules about the 
relationship between intensity distribution of patterns 
and \lm{kR)\ values. As discussed before, the high Q 
modes (low value of \lm{kR)\) show localized patterns 
near the boundary or the center of the stadium-shaped 
microdisk cavity. As \lm{kR)\ value increases, the high 
intenty spots move toward inside the cavity and emis- 
sion from cavity becomes clear, reflecting higher enegry 
loss. The modes with highest energy loss, corresponding 
to the large \lm{kR)\ values, are the bouncing ball type 
patterns. 

In usual experimental situations, it is hard to directly 
observe either the resonance pattern inside microcavity 
or the near field distribution. Therefore the far field 
distribution of resonance would be a convenient phys- 
ical quantity with which one can distinguish different 
resonance modes in experiments. Similar to the previ- 
ous pattern classification, the far field distributions of 
the resonances are also, although less obvious, classified, 
and they are roughly expected from the extension of the 
corresponding near field distributions. For example, the 
high Q modes in the odd-odd and even-odd symmetry 
classes (Fig. 3 (a),(b) and Fig. 4. (a),(b)) show good di- 
rectional emission at about 40*^ from the x— axis in the far 
field distribution, while the high Q modes in the odd-even 
and even-even symmetry classes (Fig. 5 (a),(b) and Fig. 6 
(a),(b)) have oscillatory far field distributions implying 
several emitting beams over rather wide angle range. 



IV. PATTERN CHANGE WITH REFRACTIVE 
INDEX 

In this section we study the pattern change of reso- 
nances with varying refractive index ric. Sometimes it 
seems to be accepted as a general concept that the high 
refractive index limit {ric ^ 00) in microcavity would 
be equivalent to the infinite wall case (usual billiard sys- 
tems). This might be based on the following reasoning. 
Since in higher Uc case the critical angle for total inter- 
nal reflection is very small, and approaches to zero in the 
limit, rays are confined in the cavity by the total internal 
reflections. This situation is very similar to the infinite 
wall case. One would, therefore, expect that patterns of 
eigenfunctions of the infinite wall case would represent 
the limiting {uc 00) patterns of resonances. 

Here we examine this expectation by observing changes 
of the complex wavenumber k and resonance pattern with 
increasing the refractive index Uc inside the cavity and 
by comparing the limit values with the eigenvalues of the 
infinite wall case. In order to keep the number of in- 
tensity spots in resonance patterns, we scale the size of 
the microcavity so that UcR remains constant. In Fig. 7 
is shown the evolution of resonance position kR of the 
resonance mode OO4. The Re{kR) converges to about 
4.752 and the lm{kR) seems to approach zero as ric in- 
creases. In fact, it is very difficult to say definitely if the 
lm{kR) become zero at the limit or not, due to the nu- 
merical error in BEM. The inset in Fig. 7 (b) shows the 
saturation behavior to a nonzero value of lm{kR) when 
A/" = 50, but if we increase the number of elements, the 
saturation value gets closer to the real axis with a con- 
vergence slower than logarithmic one. We also calculate 
eigenvalues of the stadium billiard in the concerned range 
which are kR = 4.23,4.44,4.81,4.86,5.12,5.46. We note 
that any of them does not coincide with the limiting value 
4.752, which mean that, although the imaginary part of 
the wavenumber kR gets closer to real axis, the resonance 
never approaches an eigenmode of the infinite wall case. 
This result can be understood in terms of the evanescent 
effect, i.e., tunneling effect, in the dielectric microcavity 
system. In the infinite wall case, the evanescent effect is 
not allowed by the Dirichlet boundary condition, and no 
amplitude outside the boundary. Although the classical 
dynamics is similar, the evanescent effect can appear in 
the microcavity case, because rays can exist outside as 
well as inside the cavity. The same argument is availible 
for the WGM in which the small emission from the cavity 
is the result of the evanescent effect. From this reasoning 
we guess that there would be very small nonzero limiting 
value of lm{kR) due to the evanescent effect. 

The resonance pattern with no /ric = 0.01, no is the re- 
fractive index outside the cavity (no = 1 in our case), is 
shown in Fig. 7 (c). The pattern inside the cavity has no 
crucial difference from the nc = 2 case when compared 
to Fig.l (a), and the emission pattern outside the cavity 
is drastically reduced, reflecting the strong ray confine- 
ment. This behavior is also true for other resonances. 



6 




(e) 



FIG. 8: (Color online) Stationary lasing modes with low 
pumping rates in the nonlinear dynamical model. The cen- 
ter positions of the gain profile are on lie(kR) = 4.3(a), 
Re{kR) = 4.7(b), Re{kR) = 4.9(c), Re{kR) = 5.1(d), and 
kR — 5.3(e). The corresponding resonances are the high Q 
modes, (a) OEi, (b) EE:^, (c) OO4 + OE^, (d) EOq, and (e) 
OO7. 



and the limiting pattern of OO5 is also shown in Fig. 7 
(d). 



V. DYNAMICAL MODELING OF MICRODISK 
LASER 

Conceptually, a dynamical description of the microdisk 
laser is the same as that of any other type of laser: it 
should include equations for material polarization and 
inversion of the active medium and equations for electro- 
magnetic field. Of these equations, only the field equa- 
tions require a special treatment in the case of the mi- 
crodisk laser. The main difficulty is that an analytical 
description of resonant mode functions, their frequencies 
and linewidths is generally impossible. Even if all these 
data are known from independent calculations, for exam- 
ple, obtained by using BEM as in previous chapters, it 
is not obvious which of these modes should be included 
into the model to incorporate the main features of the mi- 
crodisk laser dynamics. One can, of course, include into 
the model all the resonances whose frequencies fall into 
the positive gain range of the active medium. However, 
such an approach can encounter numerical problems, es- 
pecially for sufficiently large cavities and wide spectral 
contour gain. 

An alternative approach to the problem was proposed 
by Harayama et al.li]. They developed a model for the 
microdisk laser and called it " Schrodinger- Bloch model" . 
In their derivation of the Schrodinger-Bloch model equa- 
tions the authors made an approximation, replacing the 



space-dependent refractive index in some expression by a 
constant value. Using this model Harayama at al. have 
found that lasing patterns obtained are in a very good 
agreement with the BEM results. However, the authors 
did not discuss the question, why the effect of the refrac- 
tive index replacement, made to derive the Schrodinger- 
Bloch model equations, is as minor. 

Our results presented in Sec II VI allow us to clarify this 
question. Actually, the replacement mentioned above 
is equivalent, in a sense, to the refractive index change 
described in Sec II VI Fig. 7 shows clearly that the pat- 
tern and the wavelength location of the resonance do not 
change significantly with refractive index. However, the 
situation is quite different for the width of the resonance, 
i.e. for the mode damping : it is highly sensitive to re- 
fractive index variations. This behavior is quite simi- 
lar to that of conventional Fabry Perot resonator 18j, for 
which mode frequencies are determined by its length and 
mode damping - by the reflectivity of mirrors. So, we 
believe that the accuracy of the threshold characteris- 
tics obtained using the Schrodinger-Bloch model are not 
of high accuracy. However, we believe that the relative 
positions of resonances shown in Fig. 2 do not change. 

Keeping this discussion in mind, we use the 
Schrodinger-Bloch model to simulate the stadium-shaped 
microdisk laser in the oscillating regime. The mate- 
rial parameters used correspond to semiconductor active 
medium. The polarization relaxation time T2 = 50/s is 
typical for charge carrier concentration n lO^^cm"^. 
For the wavelength near A = 0.86/im and Re{kR) = 4.3 
we get the FWHM of the optical transition ARe(/ci?) = 
0.08. The transition dipole matrix element p = 15D = 
1.5 X 10~^^CGS corresponds to semiconductor quantum 
well structures. A typical value for the excited state re- 
laxation time for semiconductors is Ti = 10~^^sec, but 
we used in simulations ten or more times lower values. 
We can get various stationary lasing modes whose ba- 
sic patterns can be found in the resonance patterns pre- 
sented in Fig. 3-6. The stationary lasing modes with low 
pumping rates are shown in Fig. 8. When the gain pro- 
file is centered on Ile{kR) = 4.3, the resonance mode 
OEi is excited as the lasing mode. As shown in Fig. 2, 
it is consistent with the fact that the highest Q mode 
near Ke{kR) = 4.3 is the mode OEi. The same argu- 
ment is valid for other cases, i.e., when the position of 
the gain profile are Re{kR) = 4.7,5.1, and 5.3, the sta- 
tionary lasing modes are the high Q resonance modes, 
EEs^ EOq^ and OO7, respectively. The asymmetric las- 
ing mode shown in Fig. 8 (c) can arise from the multimode 
lasing with different symmetry classes, i.e., the resonance 
modes OO4 and OE^. The direction of asymmetric emis- 
sion oscillates, upwards and downwards alternately. We 
can expect if the pumping rate increase, this multimodes 
operation becomes single mode operation resulted from 
the locking process 15j. 

From the above calculation for the active microcavity, 
we can see that the interaction between field and active 
medium does not give any substantial pattern change 
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FIG. 9: The 50 trajectories with long path length, les > 40, 
when = 1, before escaping. 

from the resonance patterns of the passive micocavity. 
This result means that the anlysis of the resonance pat- 
terns of the passive microcavity shown in Sec lIIII would 
play an essential role in understanding the lasing modes 
from the active microcavity. 

VI. CLASSICAL CONSIDERATION BASED ON 
RAY DYNAMICS 

It is well known that the stadium system is classically 
chaotic, so there is no special structure in the portrait of 
the Poincare surface of section, showing just chaotic sea 
over the whole phase space. However, in the dielectric 
cavity case rays can escape the cavity. The path length 
before the escape should depend on the initial point in the 
phase space, and we can obtain some structure related to 
the path length in the phase space. 

From the fact that the WGM of very high Q factor 
in the circular disk are supported by the rays surviving 
infinitely through the total internal reflection, we can ex- 
pect that the high Q modes would be supported by the 
rays with long path length before the escaping by refrac- 
tion. We select 50 initial points from which the gener- 
ated ray trajectories have long path length, /gs > 40, 
when = 1, before escaping. The trajectories starting 
from the 50 initial points are shown in Fig. 9. We can see 
that the trajectories with long path length pass near the 
boundary and the center of the stadium which coincide 
with the high intensity regions in the high Q resonances. 
Similarly, the low Q resonance modes of bouncing ball 
type are supported by the rays with short path length 
before escaping. 

More generally, we can consider some distributions 



over the phase space to characterize the ray dynamical 
behavior in microcavities, e.g., the path length distri- 
bution D{s^p) and the survival probability distribution 
P{s,p,t). The path length distribution D{s^p) can be 
given by the characteristic length of the trajectory start- 
ing from an initial point Therefore, if the ini- 
tial point is on a periodic orbit {sp^pp) and the absolute 
values of associated p values are greater than the criti- 
cal value Pc of the total internal reflection, the value of 
D{sp^Pp) becomes infinity. Since the trajectoris starting 
from the initial points near the stable manifold around 
the periodic orbit would have long characteristic lengths, 
the structure of D{sp,Pp) would follow the structure of 
the stable manifold. On the other hand, when we con- 
sider an ensemble of uniform initial points over the whole 
phase space, we can trace the flow of the points in the 
phase space and after long time due to energy escape by 
emission the points with lower weight are redistributed, 
and as a result a steady probability distribution can be 
achieved. This steady probability distribution can be 
approximated by the normalized survival probability dis- 
triubtion P(s,p, t) at a long time 13]. This distribution 
would be closely related to the high Q resonance patterns 
and have the structure of unstable manifolds. 



VII. SUMMARY 

Using the BEM, we obtained all resonances in the 
range of 4 < Re{kR) < 5.5 for the stadium-shaped mi- 
crodisk, and investigated the patterns of high Q reso- 
nances for four different symmetry classes. The high 
Q resonances show high intensity distribution near the 
boundary and central part of the cavity, which is con- 
sistent with the classical ray consideration. The rela- 
tionship between Q value and resonance pattern was dis- 
cussed, and the resonances could be classified by the 
difference of the patterns even though the classical dy- 
namics is chaotic. We also confirmed that the station- 
ary lasing modes with low pumping rate in the nonlinear 
dynamical model are the high Q modes. From this res- 
onance pattern analysis, we can expect which type of 
pattern can be excited as a stationary lasing mode. 
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